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Magnetic properties of the extended periodic Anderson model 
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We study magnetic properties of the extended periodic Anderson model, which includes 
electron correlations within and between itinerant and localized bands. By combining dynam- 
ical mean-field theory with the numerical renormalization group we calculate the sublattice 
magnetization and the staggered susceptibility to determine the phase diagram in the particle- 
hole symmetric case. We find that two kinds of magnetically ordered states compete with 
the Kondo insulating state at zero temperature, which induces non-monotonic behavior in 
the temperature-dependent magnetization. It is furthermore clarified that a novel magnetic 
metallic state is stabilized at half filling by the competition between Hund's coupling and the 
hybridization. 
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Strongly correlated electron systems with degener- 
ate orbitals have attracted great interest. One of the 
popular examples is the manganite (La, Sr)Mn03.^ In 
this system, itinerant electrons in the Cg band are cou- 
pled to localized electrons in the t2g band through 
Hund's coupling; leading to a competition between an- 
tiferromagnctic (AF) correlations and the double ex- 
change ferromagnetic correlations. This yields a com- 
plex phase diagram with various types of ordered ground 
states. Other interesting examples are (Ca, Sr)2Ru04^ 
and Lan+iNinOsn+i.'^'^ In these compounds, the chem- 
ical substitution or the change in temperature is sug- 
gested to trigger an orbital-selective Mott transition,^' ^ 
where some of the orbitals become localized by electron 
correlations, while the others still remain itinerant. It 
is also proposed that in these compounds, localized and 
itinerant electrons arc hybridized with each other, induc- 
ing heavy fcrmion or bad metal behavior at low temper- 
atures.^"^" 

An important point in the above compounds is that 
the localized and itinerant bands in the d-orbital and 
their correlations play a crucial role in stabilizing the 
magnetically ordered state or heavy fermion state. Gen- 
erally speaking, in the systems with localized bands, the 
hybridization together with local electron correlations 
screen spins, leading to heavy fermion behavior with the 
large density of states (DOS) around the Fermi level, 
the so-called Kondo effect. ^^"^^ In contrast, Hund's cou- 
pling enforces parallel spins in different orbitals, enhanc- 
ing magnetic correlations, as discussed for the mangan- 
ites.^**"^^ Therefore, an interesting question arises how 
robust the nonmagnetic ground state is in systems with 
localized and itinerant bands. In a previous paper, ^ we 
have investigated the extended periodic Anderson model 
(EPAM) to clarify how the Kondo and Mott insulat- 
ing states compete with the metallic state in the para- 
magnetic case. However, a magnetic instability has not 
been discussed so far, which may be important to un- 
derstand low-tcmpcrature properties in real materials 
such as some transition metal oxides and /-electron sys- 



tems. Furthermore, the competing interactions may lead 
to nontrivial behavior in the magnetically ordered state. 
Therefore, it is highly desirable to clarify the magnetic 
properties in the system with localized and itinerant elec- 
trons. 

For this purpose, we consider a correlated electron sys- 
tem which is described by the following Hamiltonian as. 
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where cl^^{ciaa) creates (annihilates) an electron with 
spin cr(='f,J,) and band index a{=: 1,2) at the ith site, 

and Uiaa = cj^crCiaa- For the band a, tj"^ represents 
the transfer integral, V the hybridization between bands. 
The intra-band and inter-band Coulomb interactions are 
described by [/„ and U', while J denotes Hund's cou- 
pling. Finally, fj. is the chemical potential. 

To investigate the correlated electron system with one 
band localized and the other itinerant, we set the hop- 

(2) 

ping integral for the a ~ 2 band to tlj = 0, for sim- 
plicity. Then this model is regarded as the EPAM with 
not only intraband interactions but also interband ones. 
Here, to discuss magnetic properties, we make use of 
dynamical mean-field theory (DMFT).i*^-2i DMFT, 
a lattice model is mapped to an effective quantum im- 
purity, where local electronic correlations are taken into 
account exactly. The requirement that the site-diagonal 
lattice Green function is equal to that of the effective 
quantum impurity then leads to a self-consistency con- 
dition for the parameters entering the impurity problem. 
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This treatment is formally exact in infinite spatial di- 
mensions and even for three dimensions reliable results 
are obtained if non-local correlations are allowed to be 
ignored. 

When an AF instability is treated in the framework of 
DMFT,^^"^^ the self-consistency equation for the sublat- 
tice 7[= {A,B)] is represented as, 

2 



G, 



loc 7 a 



(2) 



where Qq is the non-interacting Green function for the 
effective impurity model and Gioc the local Green func- 
tion. Here, we have used the Bethe lattice with infinite 
coordination, where D is the half bandwidth for the bare 
itinerant band. Note that the self-consistency equation is 
represented only by one component of the Green func- 
tion.9'26.27 Therefore, we can introduce the effective im- 
purity model, where one of the impurity bands connects 
to the effective bath. The corresponding hybridization 

function is then defined by A-^o- = (■^)^ Gioc^a{z) 

To solve the effective impurity model quantitatively, we 
use the numerical renormalization group (NRG)^^'^^ as 
an impurity solver.'^"' '^^ This allows us to access low 
energy properties, which are particularly important in 
the Kondo insulating state. Details of the method can 
be found in literature. The hybridization function 
A^CT should be determined self-consistently through the 
DMFT condition eq. (2). 

In the half-filled system without frustration, the AF 
ground state is stabilized if intersite correlations are large 
enough. At each site, the ferromagnetic Hund's coupling 
competes with the effective AF exchange coupling in- 
duced by the Coulomb interaction together with the 
hybridization.^ Therefore, two possible spin configura- 
tions are naively expected for the AF states, which are 
schematically shown in Fig. 1 (a). A magnetization for 
each configuration may be given as m^^p^^ = mi ± m2 

where = E.(-l)''"KaT " n,ax)/(2^), P^ = 0(1) 
for i G A{B), and N is the total number of sites. Note 
that these magnetizations are not ordinary order param- 
eters characterizing the AF (I) and (II) states since they 
should be finite in both states. Nevertheless, we can dis- 
tinguish between these states: When |m^|,| — |t7t,|4^''| > 
0(< 0), the AF (I) [(II)] is realized, and a first-order 
phase transition (crossover) between the two states oc- 
curs at T = 0(r 7^ 0), where T is the temperature. Fur- 
thermore, m^^p is an important quantity since it charac- 
terizes the magnetization of ions in real materials and can 
be observed in inelastic neutron scattering experiments. 
Here, to discuss how magnetic fluctuations develop at 
low temperatures, we calculate the magnetization and 
the staggered susceptibility Xioc{= drnj^p/dhjip), where 
hAF is the staggered magnetic field. The susceptibility is 
obtained as the slope of the sublattice magnetization for 
a tiny field Haf/D ~ 0.002, which gives a good estimate 
except in the vicinity of the critical temperature. 

In this paper, we focus on the half-filled EPAM to 
discuss magnetic properties. In particular, we fix the ra- 
tio A = J/U = 0.1 and U = U' + 2Jto clarify how the 



competition between Hund's coupling and the hybridiza- 
tion affects the magnetic phase diagram at low temper- 
atures. The effects of hole doping and/or the magnetic 
field, which may yield a complex phase diagram includ- 
ing several types of magnetically ordered states, ^^'^^ are 
also interesting problems. However, these are out of the 
scope in this paper, and will be discussed elsewhere. 

In Figs. 1 (b)-(d), we show the results obtained by the 
DMFT with the NRG. When U/D = 1.0 and V/D = 0.3, 
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Fig. 1. (a) Possible spin configurations, (b), (c) and (d) show the 
staggered susceptibilities and the magnetizations as a function of 
the temperature T. Solid lines represent l?"^^! and dashed lines 
arc guide to eyes. Solid and open symbols represent the results 
for the configuration (I) and (II). 

no singularity appears in both staggered susceptibilities, 
as shown in Fig. 1 (b). This suggests that the non- 
magnetic Kondo singlet ground state is realized at low 
temperatures. We also find that magnetic fluctuations for 
the configuration (II) are enhanced at low temperatures 
and the system is close to the AF (II) state. In fact, when 
the parameters are slightly changed, the AF (II) state ap- 
pears at low temperatures. The staggered susceptibility 
diverges at a critical temperature Tjv and a spontaneous 
magnetization m^^p appears below T^r, as shown in Fig. 
1 (c). It is also found that a shoulder structure appears 
in the temperature-dependent magnetization. This im- 
plies that magnetic correlations for each band are not en- 
hanced at the same temperature (see also the inset). This 
interesting feature will be discussed later. On the other 
hand, when the Coulomb interaction and Hund's cou- 
pling are relatively large {U/D = 2.0 and V/D = 0.2), 

Xaf diverges at the critical temperature and m'^Pp ap- 
pears at lower temperatures, as shown in Fig. 1 (d). The 
AF (I) state is then realized in the ground state. 

By performing similar calculations for various model 
parameters, we end up with the low-temperature phase 
diagram at half filling, shown in Fig. 2 (a). To discuss 
the phase transitions between these states in detail, we 
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show in Fig. 2 (b) the staggered magnetization for each 
band when Y I D = 0.3. Tn the case of small 11/ D. the in- 
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Fig. 2. (a) The phase diagram of the EPAM at T/D = 4.6 X 
10""*. Solid squares, crosses and solid circles represent the Kondo 
singlet, AF (I) and AF (II) phases, respectively. Open circles 
indicate the metallic state realized in the shaded area, which 
will be discussed (Fig. 4). The phase boundaries are guide to 
eyes. Dashed line is the phase boundary obtained from the strong 
coupling limit, (b) The magnetization for each band when V/D = 
0.3. Circles and triangles represent the results at T/D = 4.6 X 
10"* and 0.015, respectively. The inset shows the results for 
V/D = 0.05. 

terband hybridization screens the local moments, and no 
magnetizations appear in both bands. Therefore, in the 
region {U/D) < {U/D)ci[^ 1.5] the paramagnetic Kondo 
phase is realized, as shown in Fig. 2 (b). The increase 
in the interaction induces magnetic moments for both 
bands with opposite signs. This implies that a continu- 
ous phase transition occurs to the AF (II) phase. Further 
increase in the interaction leads to a jump singularity in 
the temperature-dependent magnetization, at which the 
sublattice magnetizations become parallel. The first or- 
der phase transition then drives the system to the AF (I) 
state at {U/D)c2 ~ 2.1. It is also found that the jump 
singularity vanishes when the temperature is slightly in- 
creased, as shown in Fig. 2 (b). Therefore, the phase 
transition between two AF states is present at zero tem- 
perature only, while a crossover occurs at finite temper- 
atures. 

The competition between these phases may be ex- 
plained by considering the strong coupling limit [U 
oo). The system is then reduced to the Kondo neck- 
lace Heisenberg model as, H — Jinter ^ji ' ^ji + 

JintraJ2i^ii ' ^iZ, whcrC Sia = J2ss' I'^ias'^ss' Cias' , 

Jinter is the intersitc effective exchange coupling, and 
the intrasite one is represented by Jintra, instead of J. 
We note that Jinter {= /U) is always positive, while 
Jintra depends on the interactions and the hybridiza- 
tion. Its magnitude is given by the lowest singlet-triplet 
gap of the local Hamiltonian Hioc as Jintra = A£'[= 



C/(v/A2 -h4(y/C7)2 - 3A)]. In the model, three distinct 
phases appear in the phase diagram, which is schemat- 
ically shown in Ficr .3 Whpn ')Y= ^ 1^ 

AF(I) AF(II) Kondo 



Fig. 3. Ground-state phase diagram for the Kondo necklace 
model. 

the Kondo singlet phase is stabilized, where z is the 
coordination number. An increase in the intersite cou- 
pling Jinter cnhances AF correlations and a second-order 
phase transition, at last, occurs to the AF (II) phase 
at a critical value jc(= 1/2). The phase boundary ob- 
tained from the strong coupling limit, which is given 
by {V/D), = (C//4i?) V(4A + {D/U)^){%\ + {D/U)\ is 
consistent with that in the EPAM, as shown in Fig. 2 (a). 
On the other hand, the AF (I) and (II) phases are sepa- 
rated by the condition j = Q[V/D = ^/2\[U / D)], where 
the localized spins 8^2 are completely decoupled and the 
phase transition never occurs. In contrast, two bands are 
coupled through the hybridization V in the EPAM, lead- 
ing to a first-order phase transition. The corresponding 
phase boundary is in good agreement with the condition 
j = 0. When V/D = 0.05 and T/D = 4.6 x IQ-^, wc 
could not find the AF (II) phase between AF (I) and 
Kondo insulating phases, as shown in the inset of Fig. 
2 (b). This is consistent with the fact that in the weak 
coupling region the energy scale of magnetic correlations 
is fairly small and the AF state is stable only at very low 
temperatures. Therefore, we believe that in the ground- 
state magnetic phase diagram, the AF (II) phase always 
appears between the AF (I) and Kondo phases. 

Next, we discuss the finite-temperature magnetic prop- 
erties such as the shoulder structure in the magnetization 
shown before. When one concentrates on the local Hamil- 
tonian Hioci the low-lying singlet and triplet states can 
be considered to be four-fold degenerate down to a cer- 
tain temperature T* = | Ai5| . When T* is sufficiently low, 
an intersite exchange stabilizes the magnetically ordered 
state, where the itinerant band is almost fully polarized 
while nearly free localized spins appear in the other. This 
reveals that orbital-selective like features appear in the 
intermediate temperature T* < T < T/v and magnetic 
correlations in the localized band are enhanced below 
T* . Such behavior is clearly found in the vicinity of the 
phase boundary between AF (I) and (II) states, where 
T* < Tn- In fact, when U/D = 2.0 and V/D = 0.3, 
{T\Tn) - (0.03215,0.0751)) and the shoulder structure 
appears in the temperature-dependent magnetization, as 
shown in the inset of Fig. 1 (c). Furthermore, we find 
nontrivial behavior in the observable quantity |m^^|. 
Namely, when decreasing the temperature, it once van- 
ishes at a certain temperature below T*, and is induced 
again, shown as the solid line in Fig. 1 (c). This non- 
monotonic behavior may be observed in the intensity of 
the magnetic peak in the inelastic neutron scattering ex- 
periments for some transition-metal oxides with localized 
and itinerant bands such as Ca2_xSrxRuO4(0 < x < 
0.5), which then allows us to discuss the competition 
between Hund's coupling and the hybridization at each 
site. 
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Wc would like to mention another remarkable feature 
for the metal-insulator (MI) transition at half filling. In 
infinite dimensions, the Hubbard and the Kondo lattice 
model have a chance to realize the MI transition with- 
out magnetism when the system is strongly frustrated. 
However, in the system without frustration, the mag- 
netically ordered state is more stable than the param- 
agnetic metallic state and the Mott insulating state at 
zero temperature, where the charge gap always opens 
around the Fermi level. -^"^'^^ The EPAM treated here is 
not frustrated, but has competing interactions between 
the orbitals at each site, which may be referred to as or- 
bital frustration. Therefore, it is not trivial whether the 
ground state is always insulating under orbital frustra- 
tion. To clarify this, we show the DOS in Fig. 4 (a). 
When U/D < 1.5, the Kondo insulating state is real- 

Or 




Fig. 4. (a) DOS in the system with the majority spin in the itin- 
erant band when V/D = 0.3 and U/D = 1.2, 1.6,2.0,2.4,2.8 at 
T/D = 4.6 X 10"*. (b) DOS at the Fermi level as a function of 
U/D. Solid (dashed) lines represent the results for the itinerant 
(localized) band. 

ized, where the hybridization gap clearly appears in the 
localized band. In the case U/D > 2.8, the structures 
in the DOS in both bands appear far from the Fermi 
level, where the polarized ground state is realized with 
the configuration (I). On the other hand, when the sys- 
tem approaches the intermediate region, spectral weight 
is built up around the Fermi level. In fact. Fig. 4 (b) 
clearly shows that a finite DOS at each band appears 
in the intermediate region, although the a = 2 band is 
almost localized in the case (2.4 < U < 2.8). This fact 
implies that a metallic phase appears between the dif- 
ferent insulating states. Since the DOS is continuously 
changed at the metal-insulator transition points, a dras- 
tic change was not found in the magnetization shown in 
Fig. 2 (b). In our phase diagram, this magnetic metal is 
stable in the weak coupling region shown as the shaded 
area in Fig. 2 (a). We recall that when V = 0{U = 0), 
the system is reduced to the Kondo lattice (conventional 
periodic Anderson) model, where the AF (Kondo) in- 
sulator is always realized except for the decoupled limit 
U = V = Q. Therefore, we claim that the competition be- 
tween Hund's coupling and the hybridization introducing 
a kind of orbital frustration plays an important role in 
stabilizing the magnetic metal. This is in contrast to the 
other correlated systems such as the Hubbard model and 
the Kondo lattice model, where the magnetic insulating 
state is always realized, as mentioned above. 

In summery, we have investigated the extended peri- 
odic Anderson model by combining the dynamical mean- 



field theory with the numerical renormalization group. 
We have discussed the magnetic properties at low tem- 
peratures to clarify how the magnetically ordered states 
compete with the Kondo insulating state. Furthermore, 
we have for the first time found the magnetic metallic 
state in the half-filled system without lattice frustration, 
which is instead stabilized by orbital frustration. 
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